Implementation of the Non-Linear Gauge into GRACE 



G. Belanger 1 ), F. Boudjema 1 '', J.Fujimoto 2 ^ T.Ishikawa 2 '', 
T. Kaneko 3 ), K. Kato 4 ), V. Lafage 2 ), N. Nakazawa 4 ), Y. Shimizu 2 ) 

1) LAPTH, Annecy-le-Vieux F-74941, France 
; 2) KEK, Tsukuba, Ibaraki 305-0801, Japan 

3) Meiji-Gakuin University, Totsuka, Yokohama 244~0816, Japan 
4) Kogakuin University, Shinjuku, Tokyo 163-8677, Japan 

1— s 

OS 



Abstract 
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A general non-linear gauge condition is implemented into GRACE, an automated 
system for the calculation of physical processes in high-energy physics. This 
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new gauge- fixing is used as a very efficient means to check the results of large 
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scale evaluation in the standard model computed automatically We report on 
some systematic test-runs which have been performed for one-loop two-to-two 
"c^ processes to show the validity of the gauge check. 

a: 

A major part of the theoretical predictions in high-energy physics is based on pertur- 
bation theory. However the complexity increases rapidly as one moves to higher orders in 
perturbations, like when dealing with loop corrections or when dealing with many-body 
final states. In many instances calculations, if done by hand, become intractable and 
prone to error. Since perturbation theory is a well-established algorithm, it is possible to 
construct a system or software to perform these calculations automatically. 

There are several systems operating as expert systems for high-energy physics [[TJ. The 
Minami- Tateya group has developed the system named GRACE. Its structure is depicted 
in the figure. The system can, in principle, deal with the perturbative series up to any 
order. For instance all diagrams contributing to a process are generated automatically 
given the order of perturbation, and specifying the particles. However, due to the handling 
of the loop integrals, practical calculations are, for the moment, restricted to tree and 
one-loop orders. Two-loop calculation is possible only for some limited cases. The GRACE 
system can work for any type of theory, once the model file is implemented. Here, the 
model file is a database which stores all component fields and interactions contained in 
the Lagrangian of the theory. Besides the model file, peripheral parts are sometimes 
required. For instance, the structure of vertex in the theory is absent in the tool box of 
GRACE, the definition should be added. Also, soft correction factor, kinematics code, the 
interface to structure functions, parameter control section, and so forth might be supplied 
if necessary. The system is versatile enough to include new features and such added 
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components would become a part of the new version of GRACE. This report is confined 
to the calculation based on the so-called standard model. The extension to SUSY is 
presented in a separated talk.0 

In contrast to the manual computation, the theoretical prediction from an automated 
system is obtained by invoking several commands at a terminal with some elapsed time ( 
which is often rather long ), due for instance to numerical integrations over phase-space. 
Especially with an automated system it is difficult to judge the reliability of the final 
result, hence a need for a built-in automated function to confirm the results. At tree- 
level, the check of gauge invariance has been shown to be powerful. Within GRACE, the 
comparison is done between the unitary gauge and the covariant gauge and we observe 
~ 15(~ 30) digits agreement in double(quadruple) precision computation. 
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Figure 1: GRACE System Flow 



In the one-loop case, the calculation in the unitary gauge does not work well. Within 
GRACE one of the problems in this gauge is that the library containing the various loop 
integrals is designed assuming that the numerator of vector particles is g^ v . For instance, 
the library for vertex functions is implemented for the numerator of 3rd order polynomial 
in the loop momentum, while the handling of 9th order one is required in the unitary 
gauge. This not only creates very large expressions but also introduces terms with large 
superficial divergences that eventually need to be canceled precisely between many sepa- 
rate contributions. 

The non-linear gauge || is introduced to make the gauge-check possible within GRACE. 
We take a generalized non-linear gauge fixing condition for the standard model[|J: 

' (d„ - iea A„-ig cos ewPZ^W+^ + tw?- (v + 5H + iRx 3 )x^ 
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We take £w = £z = (,a = 1 so that the numerator structure of vectors is g^ u as with 
the usual 'tHooft-Feynman gauge. Any of the other parameters can then provide a gauge 
check. Then the GRACE library to compute loop amplitudes works without any change. 
When a = /3 = 5 = k = e = 0, the gauge turns to the 'tHooft-Feynman gauge. For 
instance, when a = 1 (/3 = — tan 2 ^), W^X^l (W ± x :¥ Z) vertex disappears. Though in 
the non- linear gauge we have new vertices in the ghost sector, e.g., ghost-ghost- vector- 
vector vertices, we can reduce the total number of diagrams and therefore speed up the 
calculation of many processes. Of course the Feynman rules in the model file is revised 
to take into account the modifications due to the non-linear gauge. This includes all 
the vertices at tree-level as well as all appropriate counter terms needed for any 1-loop 
calculation. As a check on the new input file, we have confirmed that the results is 
unchanged from the original standard model case when all non-linear gauge parameters 
are 0. 

We have tested the system for several tree-level processes as well as one-loop two- 
to-two processes. Here we report on the tests done at the 1-loop level. The ultraviolet 
divergence is regularized by dimensional regularization where space-time dimension is 
n = 4 — 2e. In the code generated by GRACE, the divergence is kept as a variable Cuv 
which stands for 1/s — 7£ + log47r. As a first step to check the system, we have compared, 
for a given random set of the gauge parameters, the results with letting Cuv= and 
that with Cuv ^ 0. Since the agreement is exact, the system passes the first diagnosis, 
renormalizability. Then, we proceed to check that the finite result is independent of the 
choice of set for the non-linear gauge parameters. 

Some of the processes are listed in the table below. The center-of-mass energy and the 
masses used in the computation are as follows: W = 500GeV, M z = 91.187GeV, Mw = 
80.37GeV, M H = lOOGeV, M t = 174GeV. For the regularization of the infra-red di- 
vergences, we introduce a fictitious photon mass, and its value in the computation is 
A = 10 _6 GeV. In the table, we present the value of 

! 1 oc 2 »(r /0 °p-T^ e t), 
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and LG stands for the linear gauge ('tHooft-Feynman gauge) case with Cuv=0 and NLG 
does for the non-linear gauge case with Cuv=l, a = (3 = l, 5 = k = e = 0. The latter case 
corresponds to the background gauge formalism. The number of diagrams depends on 
the choice of the gauge-fixing condition since some vertices are not present in some gauges. 
The counting of the total number of diagrams in the Table refers to all possible diagrams 
and therefore with appropriate choices of gauge this number may be reduced. The number 
of diagrams involving counter-terms insertions and that with self-energy contributions is 
denoted by CT and SE, respectively. 



Number of Diagrams da jd cos 9 [pb] 
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0.04785521591 


0.0478570236 



We note that the agreement and hence the gauge independence of the result is excel- 
lent. The accuracy is only limited by that of the evaluation of the numerical loop integrals. 
Through detailed inspections of individual diagrams we have noted that, expectedly, the 
gauge independence requires different diagrams to combine in order to produce a gauge- 
parameter independent result. In this sense, the check by gauge invariance is a powerful 
diagnostics. The comparison is still in progress for other processes at one-loop to estab- 
lish the validity of the system. We can then proceed to more complicated processes, i.e., 
one-loop two-to-three or two-to-four processes and can confirm the results of large scale 
computation by the non-linear gauge method presented here. 
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